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Abstract 

We propose a formulation to obtain the line shape of a magnetic response with dissipative effects 
that directly reflects the nature of the environment. Making use of the fact that the time evolution 
of a response function is described by the same equation as the reduced density operator, we 
formulate a full description of the complex susceptibility. We describe the dynamics using the 
equation of motion for the reduced density operator, including the term for the initial correlation 
between the system and a thermal bath. In this formalism, we treat the full description of non- 
Markovian dynamics, including the initial correlation. We present an explicit and compact formula 
up to the second order of cumulants, which can be applied in a straightforward way to multiple 
spin systems. We also take into account the frequency shift by the system-bath interaction. We 
study the dependence of the line shape on the type of interaction between the system and the 
thermal bath. We demonstrate that the present formalism is a powerful tool for investigating 
various kinds of systems, and we show how it is applied to spin systems, including those with 
up to three spins. We distinguish the contributions of the initial correlation and the frequency 
shift, and make clear the role of each contribution in the Ohmic coupling spectral function. As 
examples of applications to multispin systems, we obtain the dependence of the line shape on the 
spatial orientation in relation to the direction of the static field (Nagata-Tazuke effect), including 
the effects of the thermal environment, in a two-spin system, along with the dependence on the 
arrangement of a triangle in a three-spin system. 

PACS numbers: 67.10.Fj, 03.65.Yz, 71.70.-d, 76.20.-|-q 
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I. INTRODUCTION 



Recently, the quantum dynamics of microscopic systems have been observed due to the 
development of experimental methods. For example, the quantum mechanical magnetization 
processes of single molecular magnets (SMM) have attracted much interest. Various new 
aspects of quantum effects are seen in such systems|ll, 3, I^, Q]. 

To investigate the energy level structures of SMM molecules, electron spin resonance 
(ESR) experiments have been conducted for Mni2[5i]. The quantum tunneling effect was 
mo„,tored by a proton NMR in Fe.B Q and the dynamics of each magnetic ato,n was 
studied using NMR on Mn atoms in Mni2|8|. The temperature dependence of the ESR 
signal was also studied in Visj^. 

Complex susceptibility has been studied for a long time, and the effects of the exchange 
/or dipolar interactions between the numerous constituent spins have been clarified id . 



ana 
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12l |. In order to evaluate the ESR spectra for spatially-structured systems, theoretical 



approaches for obtaining line shapes from a microscopic view point have recently been 
proposed by focusing on the effects of the interactions between spins, using direct numerical 
evaluations of the Kubo formula along with a he,dtheo.et,ca.app.oachH 



18j. In these previous works, the line width comes from the interactions between the spins. 



which are described by the Hamiltonian system, and the line shape is given by an ensemble 
of delta-functions. The effects of contact with the thermal bath have not been studied, even 
though the thermal effect has attracted interest in studies of microscopic processes. Thus, 
an approach becomes necessary to introduce the effects of the surroundings, which cause the 
temperature dependent width of each resonant peak in the complex susceptibility. 

In order to take these effects into account, we have to study an extended system in 
contact with a thermal bath, and consider the dynamical effects from the thermal bath. For 
this purpose, the time-evolution of the reduced density operator is usually studied, which is 
obtained by projecting-out the degrees of freedom of the thermal bath. A standard formalism 



has been established for the equation of motion for the reduced density operator 



IS 
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22], which is generally called the quantum master equation. This formalism has been 



successfully applied to various fields. For example, the natural line width of a two-level (spin) 



system has been estimated 
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24| . and systems with interacting spins[25|] and nonlinear 



spin relaxation 26|] have been studied. A rapid thermal bath correlation was assumed in 
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these studies, and analyses were therefore made in the Markovian hmit. 

However, the effect of the finite correlation time of the thermal bath becomes important 
when we are interested in phenomena where the time scale of the system is comparable to 
that of the thermal bath. Then, we have to take into account the time-correlation function 
of the thermal bath and the initial condition of the density operator in the above mentioned 
master equation. In the equilibrium state of the total system, which consists of the system, 
the bath, and the interaction between them, the density operator is not given by a decou- 
pled form. Therefore, we need to take the contribution from it into account, even though 
this effect has often been ignored by assuming a factorized form of the density operator. 
In the regression theorem 27|] , we obtain the time evolution for the average of a quantity 
for the factorized initial condition and estimate the correlation function from it, which is 
good in the Markovian limit 
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30|. However, for short time phenomena, we need to 



estimate the correlation function of t he q uantity in a non-Markovian evolution, treating the 
initial correlation correctly 31, 32, 33I, 34, 35|. For example, Tanimurajs^ obtained an exact 



hierarchical formulation with a functional integral for the spectral distribution of an Ohmic 
form with a Lorentzian cutoff. Breuer and Petruccione 331] studied the effects of the initial 
correlation on the dynamics of a spin-boson system, and pointed out the importance of their 
contribution. However, they did not obtain an explicit form for the correlation function or 
the complex susceptibility as a function of the frequency. 

In the present paper, we provide a formulation for the complex susceptibility by extending 
the Nakajima-Zwanzig type of master equation without discarding the non-Markovian effect 
and the initial correlation. We derive an equation for the motion of the response function. 
Then we consider the equation of motion of the quantity [A, Weg], where Weg is the initial 
density operator for the total system and A is a system operator. We include the initial 
correlation between the relevant system and the bath, which is called the "inhomogeneous 
term" of the master equation. Since the equation is described by a time-convolution(TC) 
type of equation for the non-Markovian dynamics, the Laplace transformation can be explic- 
itly evaluated. Here we obtain a concrete form of the complex susceptibility. It should be 
noted that the obtained formula is easily evaluated, even in interacting spins, by a concrete 
numerical calculation. Moreover, by using the Hilbert-Schmidt (H-S) representation, the 
formula is compactly expressed. In the present formulation, we can include the frequency 
shift due to a system-bath interaction, which comes from the imaginary part of the memory 
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term expressed by the principal value integral of the correlation function of the thermal bath 
operators. While we present the formula up to the second order of cumulants, it could easily 
be extended to the higher orders. 

We apply the obtained formula to spin systems linearly interacting with a bosonic bath. 
For a single spin system, we study the dependence of the line shape on the type of system- 
bath coupling, e.g., the case of pure dephasing, in which only the diagonal component of 
the spin interacted with the bath, and the case of longitudinal relaxation, in which the off- 
diagonal components interacted with the bath (the non-adiabatic interaction). We find that 
the initial correlation and the frequency shift due to the memory kernel are more dominant 
in the pure dephasing case than in the non-adiabatic interaction case. Owing to the usage 
of the H-S representation, we could extend our formalism to multiple-spin systems in a 
straightforward way. For a linearly coupled spin chain, the dependence of the peak shift 
on the angle between the chain and the static field has been studied as the Nagata-Tazuke 



effect [36|]. As an example of an application to multispin systems, we study the dependence, 
including the effects of the thermal environment. We also study the relationship between 
the line shape and the geometrical configuration in a three-spin system on a triangle. 

This paper is organized as follows: We provide a general formulation of susceptibility in 
Sec. 2. The application of the obtained formula to the linear spin-boson model is given in 
Sec. 3. Discussions and concluding remarks are given in Sec. 4. 

II. FORMULATION 

In this section, we present a formulation of the complex susceptibility of a system in 
contact with a thermal bath. Generally, the linear response theory gives the complex sus- 
ceptibility in the formjl^ 

X,u{uj) = lim - / dt e-*-*-^* Tr[5^(t), AJp^g, (1) 

which describes the response of the operator to an oscillating external field conjugate to 
the operator Aj, with the frequency u. Here, fi and u are components of the operators B and 
A, respectively, and peq denotes an equilibrium state. If we consider the response in a pure 
quantum state, the time evolution of -B^(t) is given by e*^s*i?e~*^^* and peg is e~^^^/Zs, 
where Tis is the Hamiltonian of the system and is the partition function of the system 
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at a temperature T (/3 = 1/ k-sT) . On the other hand, to analyze the complex susceptibility 
under dissipation, we need to describe the time evolution of B^{t) by taking into account 
the interaction between the relevant system and a thermal bath. 

As will be shown in the next section, the dynamics in contact with a thermal bath are 
not only given by the quantum dynamics of the system, but are also affected by memory 
effects inherent in the contact with the thermal bath. The memory effect is often treated 



in the so-called Markovian approximation 



37|. This approximation is often used to study 



the time evolution of the reduced density operator of a system, which leads to the quantum 
master equation. As long as the equation has the so-called Lindblad-Kossakowski-Sudarshan 
form[37] as in the field of quantum optics, the density operator is positive definite. However, 
it has been pointed out that the Markovian approximation may violate the positivity of the 
density operator. In particular, in a spin-boson model, the breakdown of the positivity has 
been explicitly reported. A method to amend this 
kind of slippage supplement in the initial conditions 381 . 
evolution with this supplement enables us to simulate the time evolution of non-Markovian 
time evolution, but its validity is limited in a time region larger than the correlation time 



breakdown has been proposed using a 
, . Using the Markovian time 



of the thermal bath 
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Since experimental development has accelerated in recent years, we need to formulate 
a line shape theory that can correctly describe the non-Markovian effect, including the 
region of the correlation time of the thermal bath. Moreover, the term for the system-bath 
correlation at an initial time in the equation for the time evolution of the density operator 
has often been ignored. However, the importance of this term has been pointed for obtaining 



a correct description of the dynamics 
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4l|. Finally, we also need a compact formula that 



can be evaluated by a concrete numerical method. For this purpose, we will present a 
straightforward way to derive a complex susceptibility that includes the initial correlation 
as well as the non-Markovian effect. 



A. Formula of susceptibility 

We suppose that a relevant system 5* is in contact with a thermal bath R and that the 
whole system is in an equilibrium state with temperature T. Defining the density operator 
of the whole system as Weq, the linear response theory is extended to give the susceptibility 
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Xfwi^) as 

X^uiio) = lim - / dt e-^-*-^* TTs+R[B^{t), A,]Weq, (2) 

where Tis+r denotes the trace operation for the whole system. When we denote 
'Hsi'^Ri'^SR as the Hamiltonians of the systems S,R, and the system-bath interaction, 
the time evolution of an arbitrary operator for the relevant system O is determined by the 
Heisenberg equation, 

j6{t) = ^[ns + nR + HsR, d{t)] = tcd{t). (3) 

Defining the total Hamiltonian as 7i (= Ti^ + TIr + H-sr), and using the relation as 

TTs+R[B,{t),A,]W,, = TTs+R{e''''B^,A,]W,g = TTs+R.[e^^'^'B,e-^'^\A,]Weg 

= Tis+RB,{e~''''[A,,We,]}, (4) 
we can rewrite Eq. ([2]) in the form, 

X,u{uj) = lim - / cite-*-*-^*Trs+Ri?^{e-^^*[A,,iye,]} 



e- 



oo 



hm I dt e-*-*-^* TisB^paM), (5) 



with 



PAAt) = TT^e-'^'[A,,W,,]. (6) 



where Tr^ denotes the trace operation over the thermal bath. With the Fourier-Laplace 
transform f[uj] = f^dte~^'^^f(t) where f(t) is an appropriate function, we find that the 
susceptibility x^ui^^) is given by 



1 

Xf.u{uj) = -TtsB^pau[u] = Xf,ui^) - ix'^ui^)- (7) 

The above formulation shows that the procedure to obtain the complex susceptibility 
reduces to obtaining pyij,[ci;]. As shown in Appendix |X1 the time evolution of PAyif) is given 
in a form of a " master" equation by using the projection operator technique. Here, we 
define the projection operator to be P = p/jTr^. Up to the second order of the system-bath 
interaction TisRi we have 



(8) 
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where the kernel S2(t) and the inhomogeneous term \l/2(t) are given by 

S2(t - r) = TrR(-z£i)e-*^°(*-^)Q(-z£i)p^, (9) 

and 

= TrR(-z£i)e-^^°*Q[^,-po / dXHii-inX)], (10) 

respectively. In Eqs. ([9]) and ffTOj) . we used the following definitions £a;C) = j\Hk,0], and 
'^o = 'Hs+'Hr and ?ii = TY^ij . Because the whole system is assumed to be in an equilibrium 
state, we have to take into account the contribution of the initial correlation between the 
system and the thermal bath, which is represented by the inhomogeneous term \l/2(^)- 
From Eq. ([7]) and Eq. ([H]), the susceptibility Xfj,u{^^) is given by 

= ^Trs4-— — ^(pA.(0) + vl>2M), (11) 

where we define iCsO = j[Hs, O] with an arbitrary operator O. Our remaining task is to 
obtain 'E:2[uj] and \l/2[co']. For this purpose, we give concrete expressions for H2(t) and \I'2(^) 
in the next subsection. 



B. Concrete expressions for H2(t) and "^2{t) 

For simplicity, we consider the case in which the interaction between the system and the 
thermal bath is given in the form 

n, = nsR = hXY (12) 

with the system operator X and the thermal-bath operator Y. In this form, the second and 
third terms in Eq. ([8]) are given by 

/ dT'^2{t-T)pAy{T) 

Jo 

= {-'j-f f^drTvRiHu [nr{-{t - r)), e-^o(*--)p^,(r)] 

= _ r - r)XX{-{t - r)){e-^-(*--)p^,(r)} 

-$(t - r)X{e-^-(*--)p^,(r)}X(-(t - r)) 
-$(-(t - r))X(-(t - r)){e-^-(*--)p^,(r)}X 
+^{-{t - r)){e-^^-(*--)p^,(r)}X(-(t - r))X], 

(13) 
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and 

^2(t) = i--)TrB[ni,[M-t),-PA / dxx{-thx-t)] 

= i / dX^{-ih\-t){{XA^{-t)pAX{-ih\-t) - XpAX{-ih\-t)A^{-t)) 
Jo 

-{A^{-t)pAX{-ih\ - t)X - pAX{-ih\ - t)i^(-t)X)}, 

(14) 

respectively. Here, we assumed that {Y(t)) = 0, and we used Eq. (1A14I) . and definitions 
Pa = ^e"^^^ with Za = Tr^e"^^-* and 

= Tin PR Y{t) Y = {Y{t) Y). (15) 

It might be convenient to use the eigenstates of an unperturbed relevant system to obtain 
the matrix elements of pAif) in the "master" equation, Eq. ([8]). We denote the eigenstates 
of the relevant system |/) and \m) for the energy eigenvalues as Ei and Em- The (/,m) 
component of S2(t — T)pAy{j) is given by 

(/|S2(t-r)p^,(r)|m) 
= - 5^[<l>(t - r)X,,,X(-(t - r))fc,„{e-'^«(*--VA.(r)}„,„. 

k,n 

-$(t - r)X,,,.{e-^^^(*-^VA.(r)}fc,„^(-(t - r))„,^ 
-$(-(t - r))X(-(t - r)\k{e-"''^'-^^ PAAr)}k,nXn,m 
+$(-(t - r)){e-^^^(*-^)pA.(r)};,,.X(-(t - r))fc,„X„,^]. 

(16) 

We can obtain the elements for '^2{t) in a similar way. In order to evaluate the susceptibility, 
Eq. ffTTj) . we need to obtain the Fourier-Laplace transform of each element and solve the 
simultaneous equations for Eq. ([8]). We can express the equation by making use of the 
Hilbert-Schmidt (or Liouville) space, which we show in the next subsection. 

C. Transformation to Hilbert-Schmidt space 

In order to evaluate the susceptibility, Eq. (fTT]) . it is convenient to transform operators 
of the relevant system S into vectors that construct the H-S space. This is because the 

9 



Liouville operators in the Hilbert space are written as a supermatrix in the H-S space, 
which makes the evaluations easier. Defining a scalar product between operators O and Vj 
as (Vj, O) = TiVjO, the transformation from the Hilbert space to the H-S space is done by 
expanding an arbitrary operator O in the Hilbert space with a set of orthonormal operators 
Vj as 

^ = J2ny^,o), (17) 

i 

where the orthonormal condition of Vj is written as (Vj,Vj) = Sij. We can transform an 
operator in the Hilbert space to a vector in the H-S space with the set of {Vi,0). In the 
case where an operator O is written as an dimensional matrix, the corresponding vector 
in the H-S space has elements. 

Multiplication operations on a density operator in the Hilbert space are transformed to 
a supermatrix in the H-S space: When the arbitrary operators Oi and Ol are multiplied by 
the density operator p as OipOl, the product is transformed into the H-S space as 

dipdl ^Mp, (18) 

where the supermatrix jCi is symbolically expressed as 

M=Mi®M*2. (19) 

Here Aii and Ai2 correspond to matrices of the operators Oi and O2, ® denotes the 
Kornecker product, and * denotes the complex conjugate operation. When the density 
operator is written in an x A^- dimensional matrix, pAu{t) is an A^^-dimensional vector 
and A^H2(^) is an A^^ x A^^-dimensional matrix. 

Using Eqs. fllSI) and ( fT9l) . we obtain the transformation of Eq. ([8]) into the H-S space as 



J^pAuit) = —MsPAuit) + dTM^,{t-T)pAu{T)+^2{t), (20) 



which gives the susceptibility in a more straightforward way than using Eq. ([T6|) . since the 
kernel, S2(t) is written in the H-S space as a matrix 

Ms^t) = -$(t)[{XX(-t)e-i^s*} ® (e-t^s*)* - {Xe-s«s*} ® {X{-t)^ e-^^^'Y] 

+$(-t)[{X(-t)e-^^^-*} ® [X^e-^'^^'Y - e"^^^* ® {X^ X{~t)^ e-^^'^^'Y]■ 
{2l) 
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The inhomogeneous term in the H-S space is given by 
m^{t) = i / dX<^{-ihX-t) 



x{{xA^{-t)} ® {x{inx - 1)}* - {X} ® {A^{-t)x{thx - 1)}* 

® {xx{inx - 1)}* - 1 ® {xA,{-t)x{thx - t)}*)}pA, (22) 

which is an A^^-dimensional vector for the x A^-dimensional density operator. 
The z(n, m)-th component of Eq. (120|) is expressed by 

-^PA ui't)i(n,m) = —{■^■M.SPAv{t))i{n,m) 

+ / rfrA^H2('^ — ''")i(n,m),j(n',rn')P^'^('^)i("'>m') + ^2(^)i(n,m), (23) 
JO 

where z(n, m) = (n — 1)A^ + m with n, m = 1, 2, . . . A^. It should be noted that the Fourier- 
Laplace transform of the memory kernel in Eq. fl23l) is given by 

OO l-t 

dte"^"^^ / dTM-=^{t — T)i(n,m),j(n',rn')PAu{T')j{n',m') =J^-=2[^]i{n,m)J{n',rn')PAu[^]j{n';m')-, 

Jo 

(24) 

and, we have 

iujpAu[uj] - pAu{^) = -^MsPAu[uj] + Me2[^^]pAu[uj] +^2M- (25) 

Thus, we have 

PA.M =-^x(PA.(0) + ^2M), (26) 

with 

M^=[iiu + ^Ms-Ms2M-\ (27) 

which corresponds to (itu + iCs — S2[uj])^^ in Eq. ( ITTi) . All of the matrix elements of Ai^ 
are given in an explicit way, as will be shown below. The complex susceptibility in the H-S 
space is given in the form 

X^uiuj) = ^(4>-^x(PA.(0) + ^2M)). (28) 
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D. Concrete form of A^22M 

Now we obtain the matrix elements of A^H2[a;] on the basis of the eigenstates of the 
relevant system, 



-^S2 ['-^]i(n,m),j(n',m') 



N 



poo ^' 
k=l 

k=l 

with the eigenfrequency uji = Ei/h. A more explicit expression is obtained by introducing 
the spectrum of the thermal bath, as 

/oo 
dte-'^%{t). (30) 
-oo 

Using the spectrum, 

-^S2 [^]i{n,m),j{n' ,m') 

-j^ i-oo ^ 
= —— / d^' il"^^ Xn,kXk,n'Sm,m' l^pi^ + Uj' + UJkm) 
-^-"^ k=l 

+X„,„/X;^_„,Kp(u; + J + u:n'm)]J{!^') 

N 

~[y^ ^m,k^k,m'^n,n'l^p{^^ + Cj' + Unk) 
k=l 

+X„,„/X;^_„,Kp(u; + Uj' + UJnm')]Ji-^^')}, (31) 

where we define Ukm = ^k — and use the following relation, 

lim / dt e"*^*-"* = 7r5(cu) - ip- = Kp{uj). (32) 

The terms of the principal value represent the frequency shift that results from the system- 
bath interaction. While these terms have often been neglected, we can take them into 
account in the present formalism. 

E. Concrete form of 

The inhomogeneous term, Eq. ( IT^ . is simply written in the H-S space by the multiplica- 
tion of a matrix and the H-S vector of pa_ 

^2it)=M^^{t)pA, (33) 

12 



where 



■^^^2 i.^^i(n,m),j{n' ,m') 

r"/3 



2 ' 

N 



= t [ dX^-^hX-t){y2Xn,kA,k,n'X*^m'e-'^''"''e' 

k 

N 

k 

N 

— (A ,p-«'^nn'* )f* X* p*(*+*^-^Vfcm' 

k 

N 
k,l 

The Fourier-Laplace transform of the inhomogeneous term jCi^^(t) is given as 

■M.^^[u!]i(n,m),j(n',m') = {/_^ -^n,kAuk,n'Xm,m' '^pi^ + ^^kn' — ^mm') + ^mm') 

J-CXD ^ 

N 
k 

N 

— {Avn,n' Xm,kXk,m''^pi^' + ^nn' — i^km')Ki{uj' + CUfcm') 
TV 

-l^n.n' ^ Xm,kAlk,lXlm''^pi^' ~ ^kl " + ^im'))}' (35) 



where we define Kjiu) as 



/5 1 _ 

rfAe-^^'" = = Ki{uj). (36) 



The complex susceptibility, Eq. (!28|) . is now written in an explicit form, which can be applied 
to an arbitrary type of thermal bath by specifying J{uj). We will show a few examples of 
baths in the next subsection. 



Bath 



When analyzing the relaxation phenomena of a relevant system, we often introduce a 



thermal bath that consists of an infinite number of bosons 



42 



43 



44l | or spins 45|] . This 



section discusses procedures to obtain the spectra J{uj) for a bosonic bath . We use a 
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bosonic bath for the relaxation phenomena caused by phonons in a medium or photons in a 
cavity. The Hamiltonian for the boson system is written as, 

TYr = ^ huablba, (37) 

a 

where [b]^) denotes an annihilation (creation) operator for thea-th mode of a boson. As 
an example, we will consider a case in which the contribution to Tii, Eq. (1121) . from the bath 
is given by, 

Y^J^aM + ba). (38) 

a 

Then, the correlation function ( Y{t) Y ) for the bath is written as 

m = J29'M-bi)e-^--' + {bib^)e'--'}, (39) 

a 

where Qa is the coupling constant between the relevant system and the a-th mode of the 
boson. In order to evaluate the correlation function of the thermal bath, we need to introduce 
a coupling spectral function I{u) as 

Iiuj) = J29l5iuj-u;^). (40) 

a 

We can rewrite the weighted summation for an arbitrary function f{uJa) in the following 
form, 

/•oo roo 

Y.3lfM= dujJ29l5iiU-uj^)fiLu)= dujliuj)fiu). (41) 
Using Eq. fHT]) . $(t) is rewritten as 

POO 

$(t) = / du'Iiu'){{n{uj') + l)e-'"'* + n(u;')e'"'*}, (42) 
Jo 

where n{u) is the boson distribution function given by n(u;) = ^0nL_i ■ The spectrum of the 
thermal bath is obtained in the form 

/oo 
dte''^'^t) = I{uj)[n{iu) + l\e{uj) + I{uj)n{uj)e{-uj), (43) 
-oo 

where 9(u) denotes the step function. 
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III. APPLICATIONS 



We are now in a position to apply the formalism presented in the previous section to the 
relaxation phenomena in a spin system. First, we evaluate the spectral line shape of a system 
where a single spin interacts with a bosonic bath. Although such a system is trivial, the 
evaluation shows the concrete procedure, which is essentially the same as in multiple-spin 
systems. Next, we demonstrate the Nagata-Tazuke effect for two and three spin systems, 
showing the dependence of the line shape on the angle between the spatial configuration and 
the direction of the static applied field. We include the initial correlation and the frequency 
shift of the line shapes. 

A. Spin-boson model 



Suppose that a spin {S = |) linearly interacts with a thermal bath that consists of bosons. 
The Hamiltonian of the relevant system is written as 



flUJoSz, 



(44) 



and the interaction operator X in Eq. (JT2l) is given by 



X = aSx + cSz 



(45) 



where Sm-, (^^ = y, z) corresponds to the x, and z components of the relevant spin. In 
the following, we set a = sin A and c = cos A, since the generality is not lost when a is 
a real number. We control the types of relaxation by the value of A: the case of A = 
(a = 0, c = 1) describes the pure dephasing phenomena of the spin. For other cases of 
A 7^ 0, we can include the longitudinal relaxation in the transverse relaxation of the spin. 
Equation (1271) is now given as follows: The second term is written as 

/O o\ 

zcjo 
-iujo 
\0 / 

We can evaluate the third term of Eq. fl271) concretely by using Eq. fl3T|) as. 



h 



h 



(46) 
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( \a\'^(t)i[uj,ujQ] -ac(2Fs_[co',u;o] + </>3[^^,0]) -a*c(2F+[co', tuo] - (/'sftu, 0]) -|ap<^i[u;, cjq] \ 
— a*c04_ [cu, cjo] |ap044.[co', 0] + 2c^04_[(x;, tuo] — ct*^04+[t^, 0] a*c04_ [cu, c^o] 

— ac04_,_[u;, cjo] — a^04+[i^,O] |ap04^_[co', 0] + 2c^04_,_[u;, Wq] ac04^_[ci;, cjq] 

\ -|a|V2[t^,t^o]) ac(2F_[c^,u;o] - 03[t^,O]) a*c(2Fs+[u;, cjq] + 03[t^, 0]) |a| ^2!^^, t^o] / 

(47) 



where 



(PiIuJ.UJq] = F+[uJ,UJq]+ Fs^[uJ,Uq] , (j)2[uJ,UJo] = Fs+[u,Uo] + F^[u,UJo], 

(f)3[u,uo] = F+[uj,uJo] - Fs+[u,uo] , (f)4^[u,uJo] = F±[uj,uJo] + Fs±[uj,uJo] (48) 



and 



F±[uj,uJo] 



i(±LO0—Uj)t 



dt^(t)e 



7r{/(— ± ujQ){n{—uj ± ujq) + 1)6'(— c<j ± cuq) 

+I{-{-uj ± ^o))ti(-(-^ ± ^o))^(-(-^ ± ^0))} 

-zp r dio' ( ^ \ A ^')(^(^') + 1) + ^ 



U ^ Uq — LJ 



-I{uj)n{uj)} 



POO 

Jo 



(49) 



= ■k{I{uj =f u;o)(n(u; =F ujq) + l)9{uj =F cuq) 

+I{-{uj T uJo))n{-{ij T ujo))0{-{uj T ^0))} 

-zp HdJ — + 1) + ^^—-A^'X^')}) ■ 

(50) 

It should be noted that the principal value integrals are included in Eq. (!50l) . The inhomo- 
geneous term is given by 

a*cA^+r]i^[uj] + acAi,_T]i+[uj] - \a\'^A^^ + V2+[^]) ^ 

^ 2a*cA,,r]24^]+a*^A,+r]34t)-A,_{2c^r]i4uj] + \a\%4t)) 

'^2[^\ = 2% X , (51) 

2acA^,r/2-(t) + a^A^_7]3+{t) - {2c^Vi-M + |oP?73- W) 

\ -acA^_rii+[^^] - a*cA^+rii4<^] + \a\'^A^^ {V2-[^] + ^2+M) / 
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where ri^^±[u]{fi = 1, 2, 3) are given in Appendix B. 

Next, we show a numerical evaluation of the susceptibility for the Ohmic coupling spectral 
function 

/(cu) = s w 6"'^/'^% (52) 

where s denotes the coupling strength and Uc denotes the cut-off frequency. 

The concrete techniques used in the numerical evaluation of the formula, Eq. fl28|) . are 
as follows. First, we avoided the matrix inversion procedure in Eq. (1271) . By rewriting the 
complex susceptibility in the form, 

t 

Xfiui^) = ^{B^,x) (53) 

with 

x = M^{pAu{0) + ^2[uj]), (54) 

we find that the essential task required in Eq. ( l28l) is obtaining the vector x by solving the 
simultaneous equation for the elements 

[tuj + '-Ms - M^M]S = (pa uiO) + ^2M). (55) 

Second, in solving the above equation, we numerically calculated the principal value integral 
in Eqs. (1491) and (|50|) by using a Mathematica built-in function. Because this evaluation 
should be performed carefully, we checked the results by comparing them with those obtained 
by the trapezoidal numerical integration method. In the present form of I{uj), Eq. fl52|) . we 
have an analytically evaluated correlation function $(t). 

The function ^Ij'{z) = -j^i^iz) in Eq. fl56|) is defined by the digamma function, ^p{z) = 



as in 461]. Using this form, we also checked the above mentioned numerical estimations of 
the time integrals in Eqs. ( l49l) and ( l50i) . We confirmed that all three of the estimations gave 
the same result. 

Figure [1] shows the imaginary part of the transverse susceptibility x"_(cl;), found by using 
Eq. fl28l) as a function of the frequency of the external field scaled by the Larmor frequency 
of the spin uq, i.e., u = uj/uq. We scaled the coupling strength s and the cut-off frequency 
Uc with cuo, and set them to be s = 0.1 and Uc = ^^c/^^o = 0.5. We also set the temperature 
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of the bath to be fc^T = huo/b. We study the three cases of A = 0, |, |, which determine 
the types of relaxation: the case of A = corresponds to the adiabatic interaction case 
(i.e., the pure dephasing case), A = ^ to the transverse relaxation case, and A = | to the 
non-adiabatic interaction case. 
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FIG. 1: Transverse susceptibility x+-i^) for ksT = hujQ/5, s = 0.1 and uic = 0.5 with changing A 



as 0, ^, ^. 



Figure [H shows that the width of the spectra decreases with increasing A. This is explained 
as follows: for A = 0, the thermal bath affects the spin as a random magnetic field along the 
z— axis. The direction of this random magnetic field tilts toward the a;— axis as A increases 



tof. 



We can explain this fact by remembering the relaxation time obtained in the Markovian 



limit, where the transverse relaxation time T2 is given by 47|, l48l |. 

Equation fl57|) was obtained for the system-bath interaction as 

Hi = hgS-R, (58) 

where S and R are the relevant spin and bath operator, respectively. The relaxation times 
Ti and tq are given by 

^ = 2Re((/p+_ + <+), - = 2Reifi,„ (59) 

J-l To 
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where 

„2 roo roo 

¥>±T = T dte^'^''\Ri{t)R^{0)), = dt{R,{t)RM), (60) 

with R± = Rrc ± iRy. Comparing Eqs. ( fT2l) and ( l45l) with Eq. ( !58l) . we find that the case of 
A = corresponds to gR = (0, 0, gR,) = (0, 0, Y) and A = f to gR = {gR^, 0, 0) = {Y, 0, 0). 
When we consider an extreme case of {Rx(t)Rx{0)) = {Rz{t)Rz{0)) a 6{uJct), we find that 
Jr = TT- = TT— for A = 0, and = = -p— for A = ^. This means that the width of 
the transverse spectrum decreases with increasing A, which essentially explains the physical 
origin of the fact shown in FigiU 

It should be noted that the type of system-bath interaction in Eq. f l58l) is different from 
the one in Eq. f[T^ . except for cases where the thermal-bath operator R is described as 
R = Rn with an arbitrary vector n. Since the vector is written as n = z{x) in the case 
of A = (A = |), we can consider the above correspondence. While we can easily extend 
the system-bath interaction to Tii = ^j-^i^i, as in Eq. (I58l) . we chose the simple form of 
Eq. (fT2|) for a demonstration. 

In FigH], we can see the higher frequency shift due to the imaginary part of H2[ct;]. The 
detailed structures of the shifts can be seen by comparing the results with and without the 
effects of the initial correlation and the frequency shift in Figsj2] ~ |H Since the broken 
(red) lines in these figures have peaks at = 1, we find that the initial correlation and 
frequency shift cause the spectra to shift to the higher frequencies. The spectral shift and 
shape depend on the coupling strength s, and a detailed analysis of this dependence will be 
presented in a forthcoming paper. 



B. Interacting spin-boson model 

Next, we discuss the absorption spectra of interacting spins in contact with a bosonic 
bath, focusing our attention on the types of spin systems that break the Heisenberg SU(2) 
symmetry. These interactions cause shifts from the paramagnetic resonance spectra, called 
resonance shifts[36:]- Typical examples of such interactions are the anisotropic exchange 
interaction and the dipole-dipole interaction. In these cases, the Hamiltonian of the relevant 
system is as follows, 

N 

Hs = huo Si^z + Hex + Hd, (61) 

i 
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FIG. 2: (Color Online)The effects of the initial correlation and frequency shift on the transverse 
susceptibility x"_(lj) for A = 0. The other parameters are the same as in Fig(TJ The solid (black) 
line refers to the evaluation with the initial correlation and frequency shift, the dot-dash (green) line 
refers to the evaluation with just the frequency shift, the dotted (blue) line refers to the evaluation 
with just the initial correlation, and the broken (red) line refers to the evaluation without the 
initial correlation and frequency shift. In this figure, we abbreviate the initial correlation as i.e. 
and the frequency shift as f.s. 

where Si^z denotes the z component of the ith spin and Tiex is the exchange interaction 
written by 



with exchange interaction energy J, and anisotropy parameter A. Ti/j in Eq. (!6T]) is the 
dipole-dipole interaction given by 



where Vij is the vector from the spin i to the spin j, Vij = \rij\. In Eq. (!63l) . we define 
D = ^{^sY where /io is the magnetic permeabihty and yU^ is the magnitude of the magnetic 
moment that carries the relevant spin. When we consider an electron(nucIear)-spin, we have 
A'-s = gef^Bidnf^n), whcrc geidn) is the g-vsAue of an electron (nuclear) and /iB(/in) is the Bohr 
magneton (nuclear magneton). 




(62) 




(63) 
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FIG. 4: (Color Online)The effects of the initial correlation and frequency shift on the transverse 
susceptibility x+-('^) for ^ = f )• The parameters and notations of lines are the same as in Figl2j 

We can investigate the resonance shift due to these interactions by considering how the 
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spins interact linearly with a bosonic bath, as 

HsR = nXY, (64) 
^ 1 

X = ^-{a*Si^+ + QiSi-) + CiSi^z, (65) 

i 

Y ^ J^aM + ba), (66) 

a 

where we define = e*^^'' sin Ai^j and q = cos Ai_j, which control the interactions between 
the i-th spin and the bath. 



1. Two-spin system 

When the relevant system consists of two spins, the spin-spin interaction portion of 
Eq. (161]) can be rewritten as 



where 



f U U U \ ( Q \ 



hii hi2 hi3, 

^21 ^22 ^23 
hsi h^2 ^33 



'5*2,5, 



(67) 



hij 



-2(J + Do(fi^-l/3)), (^ = 1,2) 
-2{AJ + Do{Ql-l/3)) 



^ 3) 



(68) 



with 







1^^. Here we define 



ri2 
ri2 



^ sin 6i2 cos 012 ^ 



sin^^i2 sin 

COS^^12 



(69) 



where 6*12 and 0i2 are the angles of ri2 in spherical coordinates(see Fig|5]). 

We can show the dependence of the line shape on the angle ^12 when keeping 0i2 = 0. 
In FigEl we show this effect for the isotropic exchange interaction, A = 1.0, and for the 
weak dipole-dipole interaction, which is scaled by the Larmor frequency of the relevant spin 
as, Do = ^ = 0.1. The exchange interaction, which is scaled by the Larmor frequency, is 
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X 



FIG. 5: Angles in spherical coordinates for two spins 5i and 52- 

set to be J = ^ = —1. We set the scaled cut-off frequency as uJc = 0.5, and the coupling 
strength as s = 0.02. We consider the type of spin relaxation to be pure dephasing by 
setting Ai^j = A2,j = with i = 1,2. Figure M^a) shows the case of a lower temperature 
ksT = huJo/5, where we can see a sharp peak, which shows a lower frequency shift as 
9i2 increases from to | via the magic angle (= arccos(-^)). For a higher temperature, 
ksT = hujQ, we find that an additional peak appears for 612 = and ^ to give asymmetric 
spectra in Figl6]^b). In the evaluations for FiglBl we include the effects of both the initial 
correlation and frequency shift, and find the peak shifts as in the lower temperature case. 

We can explain this peak shift behavior using quantum mechanical evaluations. As typical 
examples, let us take two cases of (^^12,012) = (0,0) and (7r/2,0), for which the interaction 
Hamiltonian Eq. (1671) becomes diagonal. Defining the elements of {ha} with i = 1 ~ 3 as 
{hu, h22, h3s) = -2{Jf, Jf, Jf) respectively, we have for (^12, ^12) = (0, 0), 

Jf = Jf = J- Do/S , Jf = J + 2D0/3, (70) 
and for (^^i2,0i2) = (7r/2,0) , 

Jf = j + 2D0/3 , Jf = Jf = J- Do/3. (71) 
Using the eigenvectors of 5*^, |±), which correspond to the eigenvalues ±hujQ/2, we can 
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FIG. 6: Transverse susceptibility Xxxi^) by changing 612 from to ^ with = 0. The other 
parameters were set as lDc = 0.5, s = 1/50, Dq = 0.1, J = —1, and A = 1.0. (a) shows the lower 
temperatm'e case, k^T = hoj^/f) and (b) shows the higher temperatm'e case, k^T = Tiojq. 



obtain the eigenvalues and eigenvectors of the system Hamiltonian, Eq. ( ]6T1) . in the form, 

1 



Ea = H-Jz + K) 



\a) 



1 



Eb = h{j, - J, - jy) , \b) = —[\ + -) 



Ec = H-j, - K) 



■\{3.-h)\ + +)-{K + uJo)\--)]. 



-.[{].- Jy)\ + +) + {K + u;o)\ + +)l 



Ed = h{j, + j,+jy) , |rf) = — [I + -)-!-+)], 



(72) 



where we denote = J^^ /2 with ^ = x,y, z, and K = a/ (t^o)^ + {jx — jyY- 

In FigJTl we show the dependence of the eigenvalues Em{= Em/h) with m = a, b,c on 
the scaled dipole-dipole interaction Do for 612 = and ^12 = f- The other parameters are 
the same as in FiglHl The solid (dashed) lines refer to the energy eigenvalues for 612 = 
|). Comparing the eigenstates of the isotropic Heisenberg system, which are obtained 



'12 



in the limit of Z^o ^ 0, we can consider that the states |a), \b) and |c) correspond to the 
triplet states, |1, —1), |1,0) and |1, 1), respectively. (The state \d) corresponds to the singlet 
state, |0,0).) 
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FIG. 7: The dependence of the eigenvalues Em{= Em/K) with m = a ~ c on the scaled dipole- 
dipole interaction Dq for 6*12 = (solid line) and 612 = ^(dashed line). The other parameters are 
the same as in FigJH) 

The peaks in Figj6] reflect the transitions between the triplet states for Dq = 0.1: the 
peaks in Figj6](a) correspond to the transition between \b) and |c). The solid (dashed) arrows 
in FigHrefer to the energy differences for 6^12 = (6^12 = |) around Dq = 0.1. Due to the fact 
that the Ei, and Ec for 612 = f bend inside more than for 612 = 0, the length of the dashed 
arrow is longer than the length of the solid arrow. This means that the peak frequency for 
612 = f is lower than that for 6*12 = in Fig|6](a). We also have the other type of transition, 
between the triplet states, \a) and \b). However, the transition probability is very small in 
the lower temperature range, as in Figj6](a). The amplitude becomes larger with an increase 
in temperature, resulting in additional peaks, which correspond to the transition between 
the triplet states, \a) and \b) in Fig|6](b). 

The peak shift of ESR due to the dipole-dipole interaction in one-dimensional antifer- 



romagnets was theoretically explained by Nagata and Tazuke [36l| . They investigated the 
absorption spectra by changing the direction of the magnetic field Hq from Hq/ /c to HqJ^c. 
In experiments, the line shape is usually given as a function of Hq, but not the frequency 
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uj. Therefore, we have to evaluate the dependence of the line shape based on the strength 
of the static magnetic field, rather than the frequency, as in FiglHl The correspondence is 
discussed in the Appendix C. 



2. Three-spin system 

We now discuss the effect of dipole-dipole interaction in a relevant system with three 
spins (5* = |) that form an equilateral triangle. This is an extension of the Nagata-Tazuke 
shift to a triangle system. We can find a typical example in the antiferromagnetic triangular 



spin rings of Cu 



49| . In order to study the peak shift for these three spins, we incline the 



face of the triangle from the yz-plane to the xy-plane by increasing the angle 612 from to 
^ with = and keeping the normal of the triangle parallel to the shown in 

Figl. 




X 



FIG. 8: Arrangement of 3 spins drawn as spheres. 

Figure [9] shows the line shapes of transverse susceptibility, which include the effects of the 
initial correlation and frequency shift by system-bath interaction. Here we set the exchange 
interaction J = 1, the dipole-dipole interaction Dq = 0.1, the isotropic exchange interaction 
A = 1.0, and ksT = luq. We consider the spin relaxation type to be the pure dephasing by 
setting Ai j = A2,j = with i = 1,2. We find three peaks in the cases where ^12 = and 
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FIG. 9: Transverse susceptibility x'xxi^) when changing 612 from to ^ with (/)i2 = 0. The other 
parameters are set as lJc = 0.5, s = 1/150, Dq = 0.1, J = 1,A = 1.0, and ksT = Tiujq. 

We can explain the peak shift behavior based on the angle dependence of the energy levels 
of the relevant system. These levels consist of the lower four levels in the quartet states and 
the higher four levels in the doublet states, which are almost degenerate. Figure [10] shows 
the dependence of the lower quartet-levels (a ~ d) on the scaled dipole-dipole interaction 
Dq for 9i2 = and 612 = ^. The other parameters are the same as in FiglHl The solid 
lines (arrows) refer to the energy (differences) for 612 = 0, and the dashed lines (arrows) 
refer to the energy (differences) for 6*12 = |. These arrows are placed around Dq = 0.1, 
which corresponds to the case in Fig|9l The length of the solid arrow between |c) and \d) is 
longer than the dashed arrow between |c) and \d). Since the transition probability between 
the lower energy levels becomes higher for the relatively lower temperature, we can consider 
that the highest peaks in Fig IH] for 6*12 = and 612 = | correspond to the transition between 
|c) and \d) in the quartet states. The fact that the length of the dashed arrow between |c) 
and \d) is shorter than the length of the solid arrow between |c) and \d) shows the reason 
for the shift in the highest peaks in FigjOl Similarly, the lengths of the arrows between \a) 
and 1 6) show the shifts of the lowest peaks for 612 = and ^12 = f- 
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FIG. 10: The dependence of the eigenvalues of the quartet-levels {Em{= E^/K) with m = a ^ d 
on the scaled dipole-dipole interaction Dq. The other parameters are the same as in FiglH The 
solid line represents 9i2 = and the dashed line represents 9i2 = f • 

IV. DISCUSSION AND CONCLUDING REMARKS 

In this paper, we extend a linear response formula to include the frequency shift and initial 
correlation between the relevant system and the thermal bath. Using the projection operator 
method, we show that the time evolution of the response function can be described with a 
Nakajima-Zwanzig type of equation. We transform the equation into the Hilbert-Schmidt 
space to give a tractable formula where the density matrices are described with vectors, 
and the super operators are transformed into matrices. The obtained formula enables us 
to systematically study the line shapes for various kinds of the system-bath interactions at 
arbitrary temperatures. Moreover, with this formula it is easy to extend the relevant system 
to include multiple interacting spins. We show the line shapes for a single and for two and 
three interacting spins, which suffer from the environmental effects of a bosonic bath. 

The obtained formula enables us to evaluate the spectra, including the following three 
effects: (1) the non-Markovian effects of system-bath interaction, (2) the frequency shift by 
the system-bath interaction, and (3) the effects of the initial correlation between the relevant 
system and the thermal bath. 

While (2) and (3) are often neglected, the roles of these effects on the steady state of 
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the reduced density operator have been studied up to the second order of interaction 4^. It 
has been pointed out that these effects are essential to ensure the modification of the steady 
state by the system-bath interaction, and to prevent the steady state of the reduced density 
operator from approaching the equihbrium state of the relevant system. Since we need the 
stationary response to the external oscillating field, it is necessary to include all of these 
effects in the time evolution of the response function to obtain the line shapes. 

We apply the present formalism to spin systems interacting with a bosonic bath. For 
a single spin system, we study the dependence of the line shape on the type of system- 
bath coupling, e.g., the pure dephasing case and the case of non-adiabatic interaction. We 
find that the initial correlation and frequency shift by the system-bath interaction are more 
significant in the pure dephasing case than in the non-adiabatic interaction case. For two and 
three spin systems, we demonstrate the dependence of the line shape on the angle between 
the spatial spin arrangement and the direction of the static field. 

We compare the obtained complex susceptibility in this paper with the conventional 
one in the Born-Markovian approximation in the Appendix D. Evaluating the transverse 
susceptibility for the spin-boson model in the pure dephasing and c = 1 in 

Eq. fj45l) ). we find that the frequency shift and the initial correlation cause the considerable 
peak shiftwhich are not included in the conventional Born-Markovian approximation. 

Since the formula is written with the convolution integral as a Nakajima-Zwanzig type of 
master equation, we can systematically extend the formula to the higher orders of perturba- 
tion for the case of strong system-bath interaction 50|, ImI]- One of the authors studied this 
problem in a strongly coupled spin-boson model, and it was pointed out that the time evolu- 
tion of the transverse component of the spin was obtained in a closed form on the assumption 
of non-adiabatic system-bath interaction and a Lorentzian type of coupling spectral function, 
which enables us to systematically obtain the absorption spectra, including up to infinite 
orders of interaction. In the present formula, a similar extension is possible. Moreover, the 
transformation into the Hilbert-Schmidt space enables us to obtain a form of absorption 
spectra on demand, even for the case of the interaction of multiple spin systems in any 
spatial arrangement. We hope that our formalism will be a useful tool for analyzing the 
dynamics of various types of interacting spin systems. 
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APPENDIX A: DERIVATION OF EQ. 

We use the projection operator method to obtain the reduced dynamics of PAuif)- How- 
ever, it should be noted that the reduced density operator pAui't) is different from the ordi- 
nary reduced density operator. It is defined as p{t) = Ti^Wit) = TrRe~*^*Vr(0). Here W{t) 
denotes the density operator of the whole system, which follows the Liouville von-Neuman 
equation, 

W{t) = -iCW{t). (Al) 

In this appendix, we show that we can obtain the "master" equation for pAui^) by focusing 
on the time evolution operator e~*^* and extracting the relevant part from it. 

Let us define a projection operator V, which eliminates the variables of the thermal bath, 
to obtain the relevant part of the time evolution operator. The projection operator satisfies 
the idempotent relation, = V. We also introduce a complementary operator Q = 1 —V. 
Here we fo 
operator |19l. 
operator as 

x{t)=Ve-'^', yit) = Qe-'^\ (A2) 
with an initial time to = 0, we obtain 

j^x{t) = V{-iC)x{t) + V{-iC)y(t) , (A3) 

and 

j^yit) = Q{-zC)x{t) + Q{-iC)y{t) . (A4) 
Equation (1A4I) has the following formal solution 

y(t)= f e-Q'^'^'-^^Q{-tC)x{T)dT + e-^'^'Q. (A5) 
Jo 
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low the standard method for deriving the equation of motion of the reduced 
20, 2l|, [5^. Denoting the relevant and irrelevant parts of the time evolution 



Substituting Eq. (lASP into Eq. (lA3p . we obtain 

^x{t) = V{-iC)x{t) + V{-iC) f e-s^^(*-") Q{-iC)x{T)dT + V{-iC)e-^"^'Q. (A6) 

Jo 

We set the specific form of the projection operator to be P = p^TrR, where pr denotes the 
density operator of the thermal bath, which is in the equihbrium state. When we multiply 
Eq. (IA6I) by the initial density operator of the whole system ly(to), from the right hand 
side, we obtain the ordinary master equation for the reduced density operator p{t). Instead 
of Wito), we could also multiply Eq. flA6P by [A^,Weq\, which gives the "master" equation 
for pAui't) in the form 

^PA.(t) = -kns,PAM + / dr^t - r)pA.{r) + m, (A7) 



dt h 

where the kernel ^{t) and the inhomogeneous term are given by 



m ^ P(-2£)e-s^^*Q(-z£), (A8) 
^(t) = Vi-tC)e-^''''Q[A,,Weg]. (A9) 



(AlO) 



Using the relations 

Jo 

and VCq = CqV, the kernel ^{t) in Eq. flASp is written as 

^{t) = P(-zi:)e"*^«*QT+exp[ /" rftV^»*'Q(-z/:i)Qe"'^°*']Q(-z/:), (All) 

Jo 

which gives the expansion of the kernel ^{t) = Xli=2'^«(^)- ^'i- dSD]), denotes the 
time ordering operator from right to left. Taking up to the second order in Ci, we obtain 

f dTE2{t-r)pAAr) = f dTV{-tCi)e-''''^^'-^^Q{-zCi)pA,{r), (A12) 
Jo Jo 

where we use the relation VQ = QV = 0. 

The density operator of the whole system Weg in ip{t) is expanded as 

= ie~^(^°+^^) = |e"^^«(l - f dm^{-ih\) + ■■■), (A13) 
^ ^ Jo 

where 

= e^^«*Hie-sW«*. (AM) 
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Using Eqs. (lAlOp and Eq. (1A13|) . we obtain the expansion of ipit) = Si=i ^i(^)- The low 
order terms are given by 



^i(t) = 0, 



^2(t) = Vi-tCi)e-''^'''Q[A,,-po I dXHii-thX)] 



(A15) 
(A16) 



where po = -^e ^'^^ with Zq = Ti s+b^ To obtain Eq. (]A16p . we take up to the second 
order in Hi by expanding the partition function Z for the total system as 

Z = Tr5+Be-^(^«+«^) 

= TTs+Be~^'^'{l- dXHi{-ih\) + ■ ■ ■) 
Jo 

= Zo + Tis+Be-^^^i- / dXHii-ihX) + ■■■)• 

Jo 



(A17) 



Using Eqs. flA7p . (IA12p . and (IA16p . we obtain the "master" equation for PAv{t) in the form 
of Eq. dHD. 
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APPENDIX B: FUNCTIONS OF 7?^,±M(/i = 1,2,3) 



The explicit forms of 77^^-1- [cj] , (/i = 1, 2, 3) in Eq. (!5T|) are given as follows. 

iShuJo 7r(l - e"'^^('^=^'^o)) 

?7i,±[cj] = cosh— — { , r {I{uj±uJo){n{uj±uJo) + l)e{uj±uJo) 

Z riyLJ ± LJo ) 

+Ii-ico ± uJo))n{-{uj ± a;o))^(-(a; ± cuq))} 

-zp / rf.;- . , — A^'Knioo') + 1) - -^^^I{u')n{J)}]}, 

Jq nuj \Ui) ^ uoq — uo uo^uoq + uj J 

(Bl) 

/3B^o 7r(l - e"^^"^) 
?72,±M =e 2 { — {I{uj ±ujQ){n{uj ±uq) + l)e{uj ±uq) 

-li-iuj ± uJo))n{-{uj ± cuo))^(-(a; ± uJo))} 
-zp/ , -I{uj')in{uj') + 1) 







1 _ l3h(uj'±U]Q) 1 



h{uj' ± Wo) ± cuo + a;' 



± cuo))n(-(cu ± a;o))^(-(cc^ ± cuq))} 



00 



1 _ p-Ph{Lu'±UJo) 1 

h[(jj' ± uJo) 00 — CO ± ujq 







(B2) 



-/(-cj)n(-c<;)^(-cj)} 
-zp/ rfw-( , ^/(a;')(n(w') + l) 

— ^ iii^'M^'m, (B3) 



1 _ „/3a(w'=Fi^o) 1 

-/(wO^(^')))}- (B4) 



h{uj' =F cuo) + a;' ± cuo 

APPENDIX C: CORRESPONDENCE WITH THE ESR EXPERIMENTS 

For a two-spin system, we found that the peaks of the spectra moved to lower frequencies 
when 6*12 increased from to |. In order to compare the ESR experiment for one-dimensional 



antiferromagnets by Nagata and Tazuke 36|] , we evaluated the spectra as a function of the 



magnitude of the static magnetic field Hq for a given frequency u of the oscillating field. 
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If a peak appears at u = 'jHq + Au as a function of u (Figl6]^b)), i.e., 

UJpeak{Ho) = -fHo + Au, 

a peak in the shape of a function of u 



peak _ ^peak AiO 



(CI) 



(C2) 



7 7 

where ^''^"'^ , gives the position of the paramagnetic resonance. Therefore, the peak moves 
in the opposite direction when we give the hne shape as a function of Hq. We give an 
example in Fig{TTl where we adopted an oscillating field with a constant frequency jj^ = 2. 
As the horizontal axis of the figure, we scaled the magnitude of the static magnetic field Hq 
with the magnitude of the exchange interaction energy, Hq = Ho/'~f\J\. We set the scaled 
exchange interaction energy as J' = |j| = — 1, the scaled cut-off frequency as cj^ = yjj = 0.5, 
the scaled coupling strength as s = 0.02, and the scaled strength of the dipole interaction 
as -Dq = p| = 0.1. Since the case of 6*12 = corresponds to Hq/ /c and 6^12 = | to Hq±c, 
we found that Fig{TT] shows the same feature as the resonant shift studied by Nagata and 
Tazuke. 




FIG. 11: Transverse susceptibility Xxxi^) W changing 612 from to ^ with (pi2 = 0. The other 
parameters are set as J' = —1, io'^ = 0.5, s = 0.02, D'q = 0.1, and A = 1.0. 



APPENDIX D: BORN-MARKOVIAN APPROXIMATION 

Let us show the relation between the formula in this paper and the conventional one in the 
Born-Markovian approximation. First, in the Born approximation, after the transformation 
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of s = t — r on the memory kernel Eq. (lA12p . we replace the time evolution of pA^it ~ s) as 



f dsE2{s)pAAt-s) = f dsV{-iC^)e-'''^'^''^Q{-iC^)e'''^'pAM 
Jo Jo 



dsV[ni,[ni{-s),pAM]- 



(Dl) 



Moreover, in the Markovian limit, we assume that the correlation time of the bath is much 
shorter than that of the relevant system, which means that we make the upper bound of 
the integral in Eq. (]D1|) to be infinity [531]. The inhomogeneous term in Eq. (1A7P can be 



neglected in the Markovian limit 



30| 



dt 



In these approximations, we have the time evolution of Pav{^) the form 

ds{^{s)X[X{-s),pAM] - ^{-s)[X{-s),paMX}. 



^ PAuit) = -l-[T-(-S,pAui't)] 



which is written in the Hilbert- Schmidt space as. 



d i ^ 

-nPAuit) = -T-^S pAu{t) Markov pAu{t), 

dt h 



(D2) 



(D3) 



where Markov IS givcu by 

POO 

A^Markov = " / {XX(-s) ® 1 - X ® 

+^-s){X{-s) ® {Xt}* - 1 ® {X^X{-s)^*}]. (D4) 



The complex susceptibility in the Born-Markovian limit is obtained as 



Xf.u{uj) = -{B^, Pa 1.(0)). 



where we denote M.^ as 



My, = [iuj + -J^s - A^Markov 



-1 



(D5) 



(D6) 



For the pure dephasing case in spin-boson model, (a = 0,c = 1, by setting A = in 
Eq. (I45l) ). we obtain 

/o o\ 



M 



Markov 



</'4+[0,0] 



10 
10 
\0 / 



(D7) 
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The difference of the matrix elements of Markov with those of Eq. 0471] comes from the 
replacement of pAu{t ~ s) with e^^'^'^pAuit)- 

In the Born-Markovian approximation, the transverse susceptibility for pure dephasing 
case is given by 

-2tan\i{(3huJo) 



UJ -UJQ- f 04+ [0,0] 

From the definition of Eqs. (148!) and (1501) . we find that the part of principal value integral 
cancels in this case. This means that the frequency shift is not included in the Born- 
Markovian approximation for the pure dephasing case. 

In FigIT2l we compare the transverse susceptibility x'+-{^) for A = of the both of 
the cases Eqs. fl50|) and flD8l) . The former includes the effects of the initial correlation 
and frequency shift (solid (black) line), while the latter is given by the Born-Markovian 
approximation(double dot-dash (orange) line). We find a considerable peak shift which 
reflects effects of the initial correlation and frequency shift from the Lorentzian line shape 
in the Born-Markovian approximation. 




FIG. 12: (Color Online) Comparison of the transverse susceptibility x+-('^) for A = between 
the evaluation including the effects of the initial correlation and frequency shift (solid (black) line) , 
and the evaluation in the Markovian limit (double dot-dash (orange) line). The other parameters 
are the same as in Figl2J 
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